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We investigate the critical behavior of three-dimensional 
random-field Ising systems with both Gaufi and bimodal 
distribution of random fields and additional the three- 
dimensional diluted Ising antiferromagnet in an external field. 
These models are expected to be in the same universality 
class. We use exact ground-state calculations with an integer 
optimization algorithm and by a finite-size scaling analysis 
we calculate the critical exponents u, (5, and 7. While the 
random-field model with Gaufi distribution of random fields 
and the diluted antiferromagnet appear to be in same univer- 
sality class, the critical exponents of the random-field model 
with bimodal distribution of random fields seem to be signif- 
icantly different. 



Above two dimensions, the ferromagnetic random-field 
Ising model (RFIM) is long-range ordered for low temper- 
atures and small random fields as was proven by Imbric 
HI and also by Bricmont and Kupiainen For larger 
fields the system develops a frozen domain state || which 
has been shown to have a complex, fractal structure |Q. 
It is now widely believed that there is a second order 
phase transition from the ordered to the disordered phase 
in appropriate dimensions although in three dimensions 
a complete set of values of the critical exponents fulfilling 
the predicted set of scaling relations could still not 
be established. E. g. the value of a is still controversially 
discussed ||. 

For the replica-symmetric mean-field solution || it was 
found that the critical behavior of the RFIM depends on 
the kind of distribution of random fields. Later, also 
for random-field systems on the Bethe-lattice |l(J it was 
demonstrated that the critical behavior depends on the 
distribution of random fields. Two recent letters [ pd|Jl2| 
were addressed to the question whether this is also true in 
lower dimensions. Swift et al. [jllj found clearly different 
critical behavior for random-field systems with a Gaufi- 
distribution (G-RFIM) on the one hand and a bimodal 
distribution (B-RFIM) on the other hand in four dimen- 
sions. They could not find a clear distinction in three 
dimensions. Here, it were Angles d'Auriac and Sourlas 
p2| who found differences for the critical behavior of the 
two systems mentioned above. Especially the values of 
the correlation length exponent v they found to be sig- 
nificantly different. 

The most prominent experimental realization of a 



RFIM is often asserted to be the diluted Ising antifer- 
romagnet in a field (DAFF) (for an overview see |l3| , |l4| ) . 
This system is thought to be in the same universality 
class as the RFIM [ |15|Jl6| but if the concept of univer- 
sality is violated for random-field models the question 
arises what values the critical exponents of a DAFF have. 
Therefore, in this work we investigate these three types 
of random-field systems mentioned above in three dimen- 
sions and we determine at a time three of the critical ex- 
ponents, v, (3 and 7 numerically in order to test if there is 
possibly a violation of universality. Especially the DAFF 
is examined the first time in this way at all. 

The Hamiltonian of the RFIM in units of the exchange 
coupling constant is 



H = - y Oi<jj - 2J Bi(n 



(1) 



The first sum is over the ferromagnetic nearest-neighbor 
interactions and the spin variables <Ti are ±1. The ran- 
dom fields Bi are taken either from a GauB-probability 
distribution P{Bi) ~ exp(— (J5i/A) 2 /2) or from a bi- 
modal distribution Bi = ±A. In either distributions A 
scales the strength of the random field. 

The corresponding Hamiltonian of the DAFF is 



(*,J> 



(2) 



where we have now an antiferromagnetic nearest neigh- 
bor coupling and the ti = 0, 1 represent the dilution. 
The homogenous magnetic field A breaks the antiferro- 
magnetic long-range order and in connection with the 
dilution it acts as random field p^| , ^6[ . Note, that for 
the DAFF the value of the critical A depends on the 
dilution of the system (see also pjj for a sketch of the 
phase diagram of the DAFF). 

As was shown by renormalization group arguments 
J, the three-dimensional RFIM has a zero temper- 
ature fixed point at a finite value A c of the random-field 
width and the temperature T is an irrelevant variable. 
Hence, we can use exact ground-state calculations to in- 
vestigate the critical behavior of our systems at zero tem- 
perature. 

For our numerical investigation we used a simple cu- 
bic lattice with periodical boundary conditions and lin- 
ear lattice sizes varying from L = 10 to L = 80 for the 
RFIMs and from L = 20 to L = 120 for the DAFF. We 
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used well known algorithms from graph theory 19 2lJ 
to calculate the ground state of a system at given field 
A. The calculation works by transforming the system 
into a network [E2J, and calculating the maximum flow 
in polynomial time [ p3|p4| ] . 

This method works only for systems without bond- 
frustration, so that spin glasses cannot be treated in this 
way. For most of those systems only algorithms with ex- 
ponential time complexity are known, for example the 
Branch-and-Cut method |25|,|26|]. For readers interested 
in the field we give some additional informations: Only 
for the special case of the two-dimensional spin glass with 
periodic boundary conditions in no more than one direc- 
tion and without external field also a polynomial time 
algorithm is known |2?|] . For the general case the sim- 
plest method works by enumerating all possible states 
and has obviously an exponential time complexity. Even 
a system size of 4 3 is too large. Branch-and-Cut works by 
rewriting the problem as a linear optimization problem 
with an additional set of inequalities which must hold for 
the solution. Since not all inequalities are known a priori 
the method iteratively solves the linear problem, looks 
for inequalities which are violated, and adds them to the 
set until the solution is found. Since the number of in- 
equalities grows exponentially with the system size the 
same holds for the computation time of the algorithm. 
Here only small systems up to 8 3 are feasible . For the 
spin-glass problem approximation methods like combi- 
nations of Cluster-exact approximation ]2§| ] and genetic 
algorithms |29| , |30| are more efficient: true ground states 
pl[ up to size 14 s can be calculated Q . The basic idea 
of Cluster-exact approximation is to build sub-clusters of 
spins which exhibit no bond-frustration. For these sub- 
clusters the graph-theoretical methods used here can be 
applied, which leads to a decrease of the energy in the 
total system. Genetic algorithms work by minimizing 
many configurations of a system in parallel, keeping only 
those which have lower energies and creating new con- 
figurations by combining already existing configurations 
and flipping some spins randomly. 

We now turn back to the RFIM and the DAFF. All 
degenerate ground states of the system are given [|3| by 
a set of clusters and a binary relation defined on it. Each 
cluster is a set of antiferromagnetically (DAFF) respec- 
tively ferromagnetically (RFIM) ordered spins. These 
spins are not necessarily spatially connected. Two of 
the clusters hold the spins which have in all degenerate 
ground states always the same orientation. The relation 
describes the conditions which must hold between the ori- 



entations of the other clusters in different ground states. 
Using this description all degenerate ground states can 
be analyzed. Since all systems have a finite number of 
spins the ground state is a stepwise constant function of 
the field and, hence, this holds for the measurable quanti- 
ties as well. The steps occur whenever a cluster of spins 
flips it orientation. In |34| it was shown that for the 
DAFF more than 95% of the spins do not contribute to 
the degeneracy. For the B-RFIM even 98% of the spins 
are frozen in different ground states. This is true for all 
fields A except of the finite number of fields where the 
ground state changes, i. e. a jump in a measured quantity 
of a single system occurs. 

From the spin configurations of the ground state, we 
can calculate the magnetization m = -p- er^ (respec- 
tively staggered magnetization for the DAFF) for a given 
sample. Due to the degeneracy of the ground states men- 
tioned above the (staggered) magnetization of a certain 
system does not inevitably have a unique value. Instead, 
different degenerate ground states of a given system may 
have different magnetization values although the energy 
of the different ground states is the same, of course. Nev- 
ertheless, with our algorithm we are especially able to 
find exactly the maximum and the minimum value of m. 

In Figure [l] we show the maximum and the minimum 
absolute value of the staggered magnetization of one sin- 
gle L — 16 DAFF sample. The dilution of the DAFF 
is 50%. As discussed above, to is a stepwise constant 
function. It shows strong discontinuities (jumps) at in- 
teger values n of the field A. These jumps are due to 
the fact that all the single spins flip at A = n which are 
antiparallel to the field and, hence, have a local field of n 
generated by their n neighbors. This effect has nothing 
to do with the critical behavior but it hinders the scaling 
analysis. Therefore, for the scaling analysis of the DAFF 
we used a higher dilution such that A c is well below 1 
and additionally we used larger lattices sizes such that 
the critical region is narrow enough so that all data we 
need are for values of A < 1. 
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1 Implementation details: We used Tarjan's wave algorithm 
together with the heuristic speed-ups of Traff. In the con- 
struction of the level graph we allowed not only edges (v,w) 
with level(u>) = level (v)+l, but also all edges (v,t) where t is 
the sink. For this measure, we observed an additional speed- 
up of roughly factor 2 for the systems we calculated. 



FIG. 1. Maximum and minimum of the absolute value of 
the staggered magnetization of a DAFF (L = 16, dilution 
50%) versus field. 

Taking the average over different disorder configura- 
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tions we can calculate the order parameter M = [|m|] and 
the disconnected susceptibility Xdis = L 3 [to 2 . Here, 
the square brackets denote an average taken over up to 
180 disorder configuration for the larger system sizes and 
6400 disorder configurations for the smaller systems, also 
depending on how close the random field strength A is 
to the critical one. 



and in agreement with most of the previous work. The 
error-bars are estimated from the finite-size scaling. 




FIG. 2. Averaged maximum and minimum absolute val- 
ues of the staggered magnetization versus field for the DAFF 
(dilution 55%, different system sizes). 

For Figure |^ we calculated the average of all maxi- 
mum values as well the average of all minimum values of 
to. Although these values differ significantly we checked 
that our results do not depend on whether we take the 
maximum, the minimum or an average value of the differ- 
ent occurring (staggered) magnetization values as far as 
the scaling behavior of the order parameter is concerned. 
Hence, we decided to neglect the effect of degeneracy and 
the results presented here for the order parameter and the 
disconnected susceptibility are taken from the maximum 
values of to. Note, that the considerations above do not 
concern the G-RFIM which exhibits a degeneracy of two 
only at the jumps, elsewhere it is not degenerate at all. 

Apart from that Figure demonstrates that for a dilu- 
tion of 55% A c is well below 1 so that all data we need for 
a finite size analysis are smooth functions without steps. 

In the following analysis we use the finite size scaling 
relations 



M = L- f) / v M ((A - A c )L 1 / !y 
for the order parameter and 

Xdis = L^ X ((A - A C )L^ 



(3) 



(4) 



for the disconnected susceptibility. Figure ||| shows the 
corresponding scaling plots for the data of the G-RFIM. 

Since both quantities, M and \dis should have the 
same critical field and the same correlation length ex- 
ponent we adjusted A c and v for both scaling plots at 
the same time. From this procedure it follows A c = 
2.29 ± 0.04, v = 1.19 ± 0.08, \3 = 0.02 ± 0.01, and 
7 = 3.5 ± 0.5. These are values which are not surprising 
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FIG. 3. Scaling plot of the magnetization and the suscep- 
tibility for the G-RFIM using A c = 2.29, v = 1.19, /3 = 0.02, 
and 7 = 3.5. 

Figure [| shows the same scaling plots for the data of 
the B-RFIM. From the scaling it follows A c = 2.20±0.02, 
v = 1.67 ± 0.11, /3 = 0.0 ± 0.02, and 7 = 5.0 ± 0.4. 
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FIG. 4. Scaling plot of the magnetization and the suscep- 
tibility for the B-RFIM using A c = 2.20, v = 1.67, (3 = 0.0, 
and 7 = 5.0. 
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These values of the critical exponents differ signifi- 
cantly from those of the G-RFIM suggesting that the two 
models G-RFIM and B-RFIM are not in the same uni- 
versality class. We will discuss this later in more detail. 
First we turn to the analysis of the data of the DAFF. 
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FIG. 5. Scaling plot of the staggered magnetization and the 
corresponding susceptibility for the DAFF using A c = 0.62, 
v — 1.14, P — 0.02, and 7 = 3.4. The unsealed data are also 
shown in Figure ^. 

Performing the same analysis as before for the data 
of the DAFF results in Figure ||. Here we obtain A c = 
0.62 ± 0.03, v = 1.14 ± 0.10, (3 = 0.02 ± 0.01, and 7 = 
3.4 ± 0.4. These values are in good agreement with the 
critical exponents we found for the G-RFIM. 



1 

0.8 
0.6 
0.4 
0.2 




1 1 

L = 10 ---5K-- 
L = 20 B- 
L = 40 ---©- 
L = 80 — A- 



-2 2 
(A - A C )L^ 



FIG. 6. Scaling plot of the magnetization with the data of 
the B-RFIM but with the exponents of the G-RFIM. 

To prove even more clearly that the data of the B- 
RFIM cannot be scaled with the exponents of the DAFF 
or the G-RFIM we show in Figure g data of the B-RFIM 
which are scaled with the exponents of the G-RFIM and 



A c is best-fitted. Comparing this Figure with Figure [| 
one can see that the data collapse is clearly worse. 

The following table summarizes all results we extracted 
from our finite size scaling analysis. 





V 


(3 


7 


A c 


G-RFIM 


1.19 ±0.08 


0.02 ±0.01 


3.5 ±0.5 


2.29 ±0.04 


B-RFIM 


1.67 ±0.11 


0.0 ±0.02 


5.0 ±0.4 


2.20 ±0.02 


DAFF 


1.14 ±0.10 


0.02 ±0.01 


3.4 ±0.4 


0.62 ±0.03 



To summarize, the values we determined for the critical 
exponents v, (3 and 7 of the three-dimensional G-RFIM 
are roughly in agreement with the previous numerical 
works p5|-f37[fl] . Small deviations - as far as they exist 
- may be due to the smaller system sizes used in earlier 
numerical investigations or due to the problem of equili- 
bration of these highly disordered systems in the case of 
Monte Carlo work. 

The values for the critical exponents of the DAFF 
which we determined here for the first time within 
the framework of exact ground-state calculations agree 
within the error bars with those of the G-RFIM confirm- 
ing that DAFF and G-RFIM belong to the same univer- 
sality class [[l5 16 . Also, the value of the exponents v 
and 7 agree reasonably with experimental measurements 

Interestingly the values for the critical exponents of 
the B-RFIM deviate from those of G-RFIM and DAFF. 
This result as well as the values of v are in agreement 
with previous numerical work ]l2] , p8[ . The fact that is 
zero may suggest that the phase transition is of first order 
as it is the case for the replica-symmetric mean-field so- 
lution Q. It should be noted, however, that it was shown 
by Mezard p9| ] that there is replica symmetry breaking 
for the mean field solution of a random- field model with 
m-component spins in the limit of large m. In p2^ | it was 
also concluded from the exact ground-state calculations 
that for the three-dimensional B-RFIM the phase tran- 
sition is of first order but on the other hand real space 
renormalization yielded deviating results concerning the 
order of the phase transition (see e. g. [Q). Also, the 
value of v is even higher here. Whether the transition in 
the D-RFIM is of first order or not cannot be judged by 
our simulations. 

We should mention that our results are based on a 
finite-size scaling analysis and in principle there is the 
possibility of relevant logarithmic correction to scaling 
which possibly could also explain the deviations of the 
scaling exponents of the B-RFIM from the exponents of 
the G-RFIM and the DAFF. However, such corrections 
to scaling are expected for systems at the upper or lower 
critical dimension of a system, a case which we do not 
consider here. 

The modified hyperscaling-relation ]4l| which can be 
written in the form 7 = Dv — 2/3 where D is the spatial 
dimension (D = 3 in our case) is fulfilled by both sets of 
exponents. 
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